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abstract 

We propose a new practical method for evaluating the critical coupling constant in one- 
dimensional long-range interacting systems. We assume a finite-range scaling and define its 
exponent for the logarithm of the susceptibility. We find criticality in the form of a zeta function 
singularity. As an example, we present results for a long-range Ising model. 

Dissipative quantum mechanics has drawn attention for a long time. Recently, exper- 
iments observing quantum decoherence have enhanced the interest in such systems. We 
need a deeper theoretical understanding to describe decoherence. 

We start with a dynamical system in which the target variable is surrounded by 
many environmental degrees of freedom. pQ Then, energy flow from the target system to 
the environment might appear as energy dissipation effects. To model such behavior, we 
consider a set of infinitely many harmonic oscillators linearly coupled to the target system 
and integrate over these environmental degrees of freedom. We thereby obtain effective 
interactions for the target system variable, which are, in general, of infinitely long range 
and may effectively work as dissipation. 

Decoherence physics has been studied with a double-well potential by observing Rabi 
oscillation, which is driven by quantum coherence. Due to dissipative effects, decoherence 
appears to suppress the oscillation. This is seen as tunneling suppression, or a localiza- 
tion phase transition, due to long-range interactions. Such systems have been studied by 
many authors, typically using, canonical methods, [2] the non-perturbative renormaliza- 
tion group, [3] and sophisticated Monte Carlo simulations. [4 J It is very difficult to fully 
analyze systems of the type in question and they are often approximated by the smallest 
number of degrees of freedom, that is, a two-state approximation. In this case, the system 
is merely a one-dimensional Ising model. 

The long-range Ising model has its own long history of research. Here, we only cite 
the original works[5],[6] and one recent work. [7] If the interaction range is finite, then the 
system has no phase transition. However, when the interaction range is infinite and suffi- 
ciently strong, there can be a phase transition giving rise to spontaneous magnetization. 
Actually, it is not easy to evaluate the critical coupling constant, and this requires a large 
simulation, even in the Ising case. [8] 

The aim of this paper is to briefly present a new, practical method for evaluating the 
critical coupling constant in the case of interactions with infinitely long range. First, we 
limit the range of interactions to some finite value n, and then we solve the system precisely 
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by using an extended type of decimation renormalization group method which we call the 
block decimation renormalization group (BDRG). Then, we assume a scaling relation, 
referred to as finite range scaling (FRS). We define an exponent characterizing the range 
n dependence of physical quantities and evaluate the exponent. Using the obtained FRS 
exponent, we estimate the properties of the system in the case of an infinite range. In 
this treatment, there appears the zeta function that determines the criticality. This FRS 
method is similar to the finite-size scaling method used in simulations of finite lattice 
systems to estimate the properties of infinitely large systems. The FRS method can be 
applied to any system with infinite-range interactions in the case that the system with 
finite-range cuttoff can be solved effectively. Application of this method to double-well 
dissipative quantum mechanics will be reported elsewhere. [H] 

We consider the one-dimensional Ising model, whose action is defined by 

S = - K j a i a i+j ~ h J2 a i ' i 1 ) 

i;j>l i 

where each spin variable, <7j, takes the value 1 or — 1, and h is an external field. The 
coupling constants, Kj, assumed to be non-negative, determine the interactions between 
spins separated by a distance j. 

As a first step, we study the nearest-neighbor model, in which only K\ is non-vanishing. 
The interactions can be represented by a 2 x 2 matrix 



T (o) 



exp(i^i + h) exp(-Ki) 
exp(-Ki) exp(Ki — h) 



(2) 



We define the decimation renormalization group (DRG) transformation by integrating 
out all the even site spins and reducing effective interactions among the odd site spins. 
These effective interactions also take a nearest-neighbor form. Thus this renormaliza- 
tion procedure can be iterated, and we define the k-th. renormalized interactions T^ k ' by 
induction, 

which are interactions between spins on sites separated by a distance of 2 k . Then, we 
define the susceptibility with respect to h for the infinite system as 



(4) 



h=0 



In this simplest case, the susceptibility is exactly obtained as exp(2ifi). 

Next, we include non-nearest-neighbor interactions, but we still limit the range of the 
interaction to n, that is, Kj = for j > n. This is the first stage in the FRS method. 
The standard DRG is not effective in such a situation. To carry out the computation for 
this system, we divide spins into blocks of size n. Then, the interactions are all limited 
to "nearest-neighbor" inter-block interactions. Thus, the system is a one-dimensional 
nearest neighbor block-spin system. 

One block contains n spins, and it possesses 2 n states. The inter-block interactions 
are represented by a 2 n x 2 n matrix. A renormalization group transformation, which we 
call the block decimation renormalization group (BDRG) [2], is defined for this block-spin 
matrix. Numerical calculation of the BDRG gives a precise value of the susceptibility for 
the system with interactions of range n, which we denote x(n). 
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Here, we propose two inequalities for the logarithm of the susceptibility, 



2K [ ™ ] <log X <2K§. (5) 

Two effective coupling constants, that are appropriate for the weak and strong regions, 
respectively, are defined by 

tf^E^ K eS=Y,j K i ■ (6) 

j j 

They have been called the Oth and 1st moments of the coupling constants in references 
and have played crucial roles in determining phase structures. We conjecture that these 
inequalities hold for any set of non- negative coupling constants Kj. Though we have not 
yet succeeded in strictly proving these inequalities, we have found no exceptions in our 
calculation employing BDRG. 

Physically, we can understand these bounds as follows: The weak bound, 2K^jp , comes 
from the lst-order perturbation expansion of logx, while the strong bound comes from 
the case of an almost ordered system, for which an approximate BDRG equation provides 
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an accurate description and gives the boundary value 2K^g . [9] In case of nearest-neighbor 
interactions, both equalities hold exactly. 

Hereafter, we consider systems whose coupling constants Kj decrease with j according 
to a power low, in the form 

K 3 = - ■ (7) 

3 3 P 

Such power damping behavior is expected to appear in the case of quantum dissipation 
due to the influence of the environmental degrees of freedom. In fact, if we set p = 2, the 
resultant effective dissipation is equivalent to that in the case that there is a friction force 
proportional to the velocity. [I] For this reason, this case is called the Ohmic case. For the 
above coupling constants in Eq. ([7]), the inequalities read 

277C(p)<logx<2^C(p-l) , (8) 

where ( is the standard zeta function. For example, the result of the BDRG calculation 
of the susceptibility for n = 11, p = 1.8 and r\ = [0, 1] is displayed in Fig. 1 with the lower 
and upper bounds defined in Eq. (jH]). That figure clearly shows that the logarithm of 
the susceptibility is accurately described by these boundary lines in the weak and strong 
coupling regions, respectively, and it moves from the lower bound to the upper bound 
rather quickly. 

Now, we describe the FRS method. We calculate the change in the logarithm of the 
susceptibility when the range of the interaction changes by 1: 

A(n,p,T)) = — [logx(n) -logx(n- 1)] . (9) 

At the boundaries, this should behave as 

aw = (i)' , A is] = (ir . (10) 
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Figure 1: The susceptibility for p = 1.8, rj = 

m ii j „ j i Figure 2: The FRS exponent /3 for p — 1.8, 

[0, 1 and n = 11, with the upper and lower ° , ^ ^ 

bounds. ?7 = [0, lj and n = 3, 5, 7, 9, 11, oo. 



Considering these boundary forms, we define the FRS exponent for A by 

^ M = ita !2I^2<> . 



n^oo 



We have not proved the existence of this limit. It is plausible, however, that it does exist. 
Then, the infinite-range contribution to the logarithm of the susceptibility is represented 
by the zeta function, 

Km logx(n) oc (((3(p,r})) • ( 12 ) 

As a real function, the zeta function ((z) is finite for z > 1 and diverges as z — > 1 + . Thus, 
(3{p,rj) = 1 corresponds to criticality. 

On the other hand, from the inequalities (jHJ), we have the inequalities for the FRS 
exponent 

P>i3{p,r])>p-1. (13) 

Because x(n) is monotonic with respect to n[TD], /3{p,T]) is monotonic with respect to r], 
and /3(p, rj) changes from p to p — 1 as rj changes from to oo. 

Taking account of these properties of the FRS exponent, we can draw some conclusions 
about the properties of the phase transition. For p < 1, the susceptibility cannot be finite 
as n —>■ oo for any rj ^ 0. Thus, the system is always ordered for any rj ^ 0. By contrast, 
for p > 2, the susceptibility is always finite as n —>■ oo for any finite rj. Thus, in this 
case, the ordered phase never appears. For the intermediate region satisfying 1 < p < 2, 
there can be a phase transition point at which j3(p, rj) crosses the value 1. These basic 
properties have been known for many years, and they were originally demonstrated by 
other reasoning or more rigorous arguments. [5] [6] Our aim here is to evaluate the critical 
coupling constant r] c (p). 

Now we report the results of numerical calculations of the FRS exponent to determine 
rj c . In the actual evaluation of the exponent, we employed the formula 

logA(n,p,77) -logA(n - \,p,rj) , . 

(3{n,p,r)) = — (14) 

logn — log(n — 1) 

to define f3(n). Analyzing its n dependence, we estimate the value of (3 in the n — > oo 
limit by linearly extrapolating with respect to 1/n. With the computing resource of a 
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desktop computer, the size n = 11 is the practical limit for an overnight calculation, and 
hence, in this paper, we consider data up to n — 11. 

In Fig. 2, we plot an example of the n dependence analysis for (3{n,p = 1.8, rj), with 
n = 3,5, 7, 9, 11, oo (which is the extrapolated value from n — 10 and 11). According to 
the FRS exponent inequalities, (3 changes from the weak limit, 1.8, to the strong limit, 
0.8. In the region of this monotonic decrease, there exists the point (3=1, corresponding 
to the critical coupling constant 7] c , which is found to be about 0.41. From the data, the 
decrease is not strictly monotonic, and the lower bound is broken, especially when we 
consider linearly extrapolated values. However, we believe that this is spurious behavior, 
due to the fact that n is too small. 

We should also note that around the most important region, (3 ~ 1, (3{n) is nearly 
independent of n, which is rather surprising, since the critical point 77 c is correctly found 
even with n = 3. This stability near the criticality is seen for all values of p. Contrastingly, 
linear extrapolation is necessary for the weak and strong regions, and this makes the 
transition sharper. We conjecture that at the critical point, the function f3(rf) is singular 
with an infinite derivative and may in fact jump to the strong limit value. 

Calculating the FRS exponent (3 for n = 9,p = [0.9, 2.1], 77 = [0,1], we obtain the 
contour plot given in Fig. 3. We plot (3{n = 9) without extrapolation. The thick curve 
represents the contour of (3 = 1, and the adjacent curves represent the contours with 
spacings of 0.2. The (3 = 1 curve is the phase boundary in the p-rj plane. The lower side is 
the symmetric phase, and the upper side is the symmetry broken phase. We see that the 
boundary values of the (3 inequalities accurately approximate the weak and strong region 
behavior for all p. We also see that non-monotonic behavior is stronger for smaller p, but 
it does not seem to affect the critical region. 




Figure 3: Contour map of the FRS exponent Figure 4: The critical coupling constant rj c 
[3 for n = 9,p = [0.9, 2.1], 77 = [0, 1]. as a function of p. 

Values for the critical coupling constant, r] c (p), are listed in Table I. For p = 2, 
it is not easy to evaluate t] c with high precision, and therefore here we list results for 
p < 1.99. Figure 4 compares our results with the strict lower bound obtained by Dyson 
and Griffiths [6| fTT](DG). t] c (p) > l/(2((p)), and with recent high precision Monte Carlo 
simulation results (MC) for p = 2, i] c = 0.6551(6). [8] It is seen that our results do not 
go beyond the DG bound, and actually, they are very near to each other in the region 
of small p. The MC results and our results appear to be consistent, but the region near 
p = 2 should be examined in more detail. 
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Table 1: The critical coupling constant f] c {p) 
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Vc 


V 


Vc 


V 


Vc 


V 


Vc 


1.10 


0.0474(1) 


1.60 


0.281(1) 


1.91 


0.512(1) 


1.96 


0.582(1) 


1.20 


0.091(1) 


1.70 


0.340(1) 


1.92 


0.524(1) 


1.97 


0.602(1) 


1.30 


0.135(1) 


1.80 


0.410(1) 


1.93 


0.537(1) 


1.98 


0.626(1) 


1.40 


0.180(1) 


1.90 


0.501(1) 


1.94 


0.550(1) 


1.99 


0.657(2) 


1.50 


0.228(1) 






1.95 


0.565(1) 







In this paper, we do not compute the critical exponents of the phase transition. The 
((z) function has a simple pole at z — 1, and we can deduce the susceptibility critical ex- 
ponent and its transition type (standard or Kosterlitz-Thouless) by examining the details 
of the /3(p,r)) behavior near i] c . 
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